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Abstract 

We present a combined analytical approach and numerical study on the stability of a ring bound 
to an annular elastic substrate, which contains a circular cavity. The system is loaded by depres¬ 
surizing the inner cavity. The ring is modeled as an Euler-Bernoulli beam and its equilibrium 
equations are derived from the mechanical energy which takes into account both stretching and 
bending contributions. The curvature of the substrate is considered explicitly to model the work 
done by its reaction force on the ring. We distinguish two different instabilities: periodic wrinkling 
of the ring or global buckling of the structure. Our model provides an expression for the critical 
pressure, as well as a phase diagram that rationalizes the transition between instability modes. 
Towards assessing the role of curvature, we compare our results for the critical stress and the wrin¬ 
kling wavelength to their planar counterparts. We show that the critical stress is insensitive to the 
curvature of the substrate, while the wavelength is only affected due to the permissible discrete 
values of the azimuthal wavenumber imposed by the geometry of the problem. Throughout, we 
contrast our analytical predictions against finite element simulations. 
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I. INTRODUCTION 


Wrinkling is a stress-driven mechanical instability that occurs when a stiff and slender 
surface layer, bonded to a compliant substrate, is subject to compression. This universal in¬ 
stability phenomenon is found in numerous natural and technological/engineering examples, 
over a wide range of length scales, including: carbon nanotubes [I], pre-stretched elastomers 
used in flexible electronics applications [2), human skin [3j, drying fruit [3], surface morphol¬ 
ogy of the brain [5] and mountain topographies generated due to tectonic stresses mm- 

Over the past decade, there has been an upsurge of interest in the study of the mechan¬ 
ics of wrinkling, along with a change of paradigm in regarding surface instabilities as an 
opportunity for functionality, instead of a first step in the route to structural failure mm- 
The first mechanical studies of wrinkling were motivated by the stability of sandwich panels 
ms. used in lightweight structural applications, in which the core acts as a soft substrate 
for the much stiffer skin. More recently, HP showed how the wrinkling of a thin film on 
an elastomeric substrate can be used to produce complex self-organized patterns. Their 
seminal work has instigated the realization of wrinkling through several different actuation 
mechanisms, including thermal mismatch na, tissue growth/atrophy [5lH3l 13] . swelling by 
a liquid pj)] or vapor solvent [16] , and pneumatics D3|- The opportunities in applications 
opened by such a wide range of external stimuli have enabled the usage of wrinkling in 
photonics ]18j . optics [15] . self-assembly [SO], microfluidics [21] and morphogenesis [22]. 

In order to provide a theoretical background to these recent developments, several authors 
have built on the pioneering work of ms. who first provided close form solutions for the 
critical stress and wavelength obtained when an initially straight beam, adhered to an infinite 
plane substrate, is placed under a state of uniaxial compression. [23] extended this work to 
consider the case of a plate adhered to a flat substrate under equi-biaxial compression and 
performed a nonlinear analysis of the Fdppl-von Karrnan equations [23!, 25]. [26] further 
refined these efforts by considering the effect of a finite substrate. Both studies showed 
the existence of multiple buckling modes associated with the same value of critical stress. 
The stability of these modes under different loadings conditions has been addressed by [27] 
[29], who produced a stability diagram covering the evolution from low to high values of 
overstress. However, experiments by [30] found disagreement at low values of overstress, 
suggesting that a finite intrinsic curvature of their experimental system, even if small, may 
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play an important role in dictating pattern selection. 

Early studies of wrinkling on curved substrates, as in the flat configuration, were also 
motivated by a structural problem; in this case, in the context of the stability of the outer 
shell of rockets [3T1 STTTTj . More recent studies that consider instabilities as a possible source of 
functionality have led to applications of curved configurations in adhesion [33j , microfluidics 
[35] . morphogenesis of microparticles [36] . optics [37] and aerodynamic drag reduction H7J. 
Curvature also plays a relevant role in the growth of biological systems [3E] • Despite these 
important emerging applications, the mechanics of wrinkling on curved substrates remains 
poorly understood, when compared to the planar counterpart. 

Systematic Finite Element simulations of wrinkling in curved systems have been per¬ 
formed MEM that highlighted a complex pattern formation process. These numerical 
studies also suggested the possibility for curvature to affect the selected patterns and modify 
the relevant characteristic length scales, which calls for a robust theoretical backing. Ana¬ 
lytical predictions are challenged by the difficulty of modeling the stiffness of the substrate, 
even in two-dimensional configurations. [CEj and [3HJ used the stiffness provided by |10j for 
the flat case, such that their model therefore neglects the contribution of curvature on the 
response of the substrate. 0] used the prediction provided by [32], which accounts for cur¬ 
vature but does not consider its influence on the wrinkling wavelength and their prediction 
does not converge to the classical planar case when the curvature tends to zero. As such, 
there is a need to quantify the effect of curvature on the stiffness of the substrate and its 
subsequent influence on wrinkling. 

Here, to the best of our knowledge, we provide the first analytical work that accounts for 
both the curvature of a (2D) shell-substrate system, as well as the finite size of the substrate. 
As an initial step, we focus our study on a curved film adhered to a cylindrical substrate, 
instead of dealing with non-zero Gaussian curvature geometries, which is left for a future 
study. We assume axial-symmetry to further simplify the system to the 2D problem of a ring 
on an annular substrate. Mechanical loading is applied by depressurizing a circular cavity 
inside the substrate, which places the system under a state of compression. This geometry is 
motivated by recent experiments ca that demonstrated the usage of wrinkling on spherical 
samples for switchable and tunable aerodynamic drag reduction. In our simplified 2D system, 
we solve the elasticity problem for the substrate and derive a close form expression for its 
stiffness, which is then used in the stability analysis of the ring to quantify the buckling 


3 


patterns. 


The paper is organized as follows: In § [TTJ we introduce our system along with its material 
and geometrical parameters. We also describe the possible instability modes, and present a 


simplified phase diagram, with the aim of providing physical insight on the problem. In § III 


we then introduce the kinematics of the ring attached to the substrate and determine the 
stiffness of the substrate. We proceed by defining a strain energy that includes both bending 
and stretching of the ring, as well as the effect of the substrate. Energy minimization yields 
the equilibrium equations of the problem. An asymptotic expansion is then used to calculate 


the principal solution and the bifurcation at the onset of instability. In § IV we describe 
the finite element simulations that we have performed for this same system. 

The results of our investigation are presented in §[VJ Throughout, we directly compare the 
analytical predictions to the numerical simulations. We start with the fundamental solution 
and the critical conditions that lead to instability. We then construct a phase diagram 
which rationalizes the dependence of the instability modes on the governing parameters. 
The results for our system are then quantitatively compared to those for wrinkling of a film 
on a planar substrate, highlighting the effect of curvature. Finally, § |VI| summarizes our 
findings and provides perspectives for potential extensions of our work in future studies. 


II. DEFINITION OF THE PROBLEM 

We study the stability of a thin elastic ring, bound to an equally curved 2D substrate 
that contains an inner cavity, a schematic diagram of which is presented in Fig. 0 a). The 
system is initially at equilibrium, with identical pressures inside and outside of the sample. 
Motivated by recent experiments on spherical specimens C2, the system is then loaded 
by applying a depressurization, P, to the inner cavity. The thickness of the ring is H , 
its Young’s modulus E F and its Poisson’s ratio v F . We refer to E F = E F /( 1 — u F 2 ) as the 
reduced Young’s modulus of the film. The substrate is made of a linearly elastic material with 
Young’s modulus E s , Poisson’s ratio v s and reduced Young’s modulus E s = E s /( 1 — I's 2 ). 
The thickness of the substrate is R — Rq, where Ro is the radius of the inner cavity. 

For convenience, we now introduce new rescaled quantities to reduce the number of 
parameters of the problem. As such, we use R and E F to normalize lengths and pressures 
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FIG. 1. (a) Schematic diagram of our system: a ring is bound to a curved substrate which contains 
a circular cavity. The system is loaded by applying a pressure differential between the inside of the 
cavity and the outside of the ring, (b-d) Representative examples of the three possible instability 
modes of a ring on a curved substrate which contains a cavity that is depressurized, (b) Wrinkling 
mode (h = 10 -2 and £ = 10 3 ), (c) global buckling mode (h = 10~ 2 and £ = 10 6 ) and (d) Biot 
mode (h = ICR 2 and £ = 10 2 ). (e) is a zoom in of (d) that exhibits the deformation of the surface 
of the inner cavity in the Biot mode. The adjacent colorbar applies to pictures b-e) and refers to 
the maximum principal component of the strain tensor of the mode, which has been normalized 
by the maximum value of each configuration. 


and define 


/, — — g — c — , — JL 

~~ R’ P ~ R ’ ^ _ E5’ P ~ Kp' 


( 1 ) 


as the dimensionless thickness, cavity size, stiffness ratio, and pressure, respectively. 

The principal solution corresponds to an axisymmetric deformation that leads to a de¬ 
crease of both the inner and outer radii. As the depressurization increases, the onset of 
instability is reached. In Fig. [ljb-e) we show representative results obtained from Finite El¬ 
ement Modeling (FEM), of the three possible instability configurations of the ring-substrate 
system, for different values of the dimensionless ring thickness, h, and ratio of stiffness, 
£. For these results, all the other mechanical properties were kept constant: the Poisson’s 
ratios of the film and substrate are vp = vg = 0.5, and the dimensionless size of the cavity 
is (3 = 0.2. In the simulations, this is achieved by fixing R = 100 units of length and Eg = 1 
units of pressure, while changing the values of H, R 0 and Ep accordingly. More details of 


our numerical simulations are provided in IIV 
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The first mode, a representative example of which is shown in Figure [jjb) for h = 10 -2 
and £ = 10 3 , corresponds to periodic wrinkling of the film with a well defined wavelength. 
The displacements are localized in the region close to the film. The second instability mode, 
for example at h = 1CU 2 and £ = 10 6 in Figure [jjc), corresponds to a global buckling of the 
structure, where both the ring and the cavity deform into an ellipse such that the wavelength 
is A = 7 tR. 

In addition to these two instability modes (wrinkling and global buckling), we have also 
numerically observed an instability on the inner surface of the cavity. However, this third 
mode does not affect the ring and is only found for low values of cavity size and stiffness 
ratio (£ = 10 2 and h = 10~ 2 in Fig. [ljc-d) . This instability was first discussed by [33] and 
we therefore refer to it as the Biot mode ; it is local in nature and only depends on the 
compressive strain at the inner surface. This type of instability mode has been recently 
studied in the case of elastomeric materials with voids by [33] and [33]. Understanding the 
specifics of this Biot mode is however outside the scope of our work and we shall not take it 
into account in our analytical model and systematic numerical investigation. 

A schematic phase diagram of our system is provided in Fig. [2j For low values of h and 
£, the ring wrinkles with a short wavelength. As either h or £ are increased, the wavelength 
also increases. Once these parameters reach a critical value, represented by the dashed line 
in Fig. [2j the instability transitions from wrinkling to global buckling. In what follows, we 
focus on rationalizing how the wavelength of the wrinkling mode, and the threshold value 
for the transition to global buckling, evolve with the elastic and geometrical parameters of 
the system. 


III. ANALYTICAL MODEL 

The ring is treated as an Euler-Bernoulli beam. The effect of the substrate is modeled as 
a restoring force that acts on the ring and is determined by solving the elasticity problem 
of the substrate with adequate boundary conditions. Minimization of the potential energy 
provides the equilibrium equations of the problem, which are solved using an asymptotic 
expansion that yields the principal and bifurcated solutions. 


6 



I tt - 

II 


vjlP 

O' 

+-> 

2 

to 

to 

<L» 

C 

tfc 

1/1 




Wrinkling 





- 1 - 

2.5x10 


Global buckling Normalized 



Dimensionless thickness, h=H/R 


FIG. 2. Schematic phase diagram of the instability modes of the system, obtained for a critical 
value of the pressure differential: i) wrinkling of the ring or ii) global buckling of the structure 
(shaded regions). The primary parameters that govern this transition are the stiffness and thickness 
ratios: £ = Ep/Eg and h = H/R, respectively. The depicted examples from FEM simulations are 
for an incompressible film and substrate, vp = vg = 0-5- They were obtained for R = 100 units 
of length and Eg = 1 units of pressure, while varying H and Ep. The colorbar refers to the 
maximum principal component of the strain tensor of the mode, which has been normalized by the 
maximum value of each configuration. 

A. Kinematics, energy formulation and equations of equilibrium 

We model the ring as an extensible Euler-Bernoulli beam made of an homogeneous and 
isotropic material. Polar coordinates are used to track the position of the ring center-line, C. 
The initial configuration of the ring, prior to depressurization, is assumed to be circular. The 
origin, O, is located at the center of the cavity, and the initial and equilibrium configurations 
of an arbitrary point of C are represented by M 0 and M, respectively. Vectors are expressed 
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in the physical base (e r , e#), derived from the polar coordinates ( r,9 ). The initial position 
of C is OM 0 = (R, 0), as shown in the inset of Fig{l|a). When the system is loaded by 
depressurizing the cavity, C deforms into a new configuration given by the position vector 
OM = R (1 + v (d) ,u (&)), where v and u are the dimensionless radial and orthoradial 
displacements, respectively. 

The infinitesimal arclcngth of C in the initial and deformed configurations are denoted 
by dso = |dOM 0 | and ds = |dOM|, respectively. Moreover, defining the tangent vector 
T = dOM/ds, allows us to express the curvature of C in the deformed configuration as 
k/R= |dT/ds|. Here, k is dimensionless and can be written in terms of v and u as 

k = 1 + (—1 + 2 u! + 2v) v" + v 2 — v — u' 2 + (—u + v') u" + ^ ^ v' 2 — u 2 ^ + li.o.t., ( 2 ) 


where the prime notation represents derivation with respect to 6 and high order terms 
(h.o.t.) are neglected under the assumption of small displacements and moderate rotations. 
We now define the elongation of the ring as e = ds/dso to express the stretching deformation, 
77 = (e 2 — l)/ 2 , in terms of v and u as 


n 




+ (1 + v) v! + V, 


( 3 ) 


so that the hoop stress in the film is Co = Eprj. 

Following Euler-Bernoulli beam theory [25], the total energy of deformation £ of the ring 
is the sum of a stretching energy £ s and a bending energy S B , 


2 ? tR 

£ = £s + £b = [ £dso, 


with 


£5 = 

£r = 


2nR_ 


E f H 2 


V ds 0 , 


0 

2ttR_ 


EfH \ 1^ 
~2AR?^ K ~ 1 dS °’ 


( 4 ) 


(5a) 

(5b) 


and £ is the energy of deformation per unit length of the initial configuration of the ring. 

2tcR 


Assuming that the reaction force of the substrate derives from a potential f Wdso, the 






equilibrium states of the ring are the solutions of 


2n R 

8£ s + 5£ b - J 8Wds 0 = 0, (6) 

o 

where 8 A is the variation of quantity A, for an arbitrary displacement held R (8v, 8u ), which 
is 27 r periodic. The computation of the variations in Eq. (J6]) leads to the Euler-Lagrange 
equations for the equilibrium of the ring, 


d£ fd£\ (d£\ _88W 

dv \ dv' J \dv” J 88v 

d£ (8£ V / 8£ \ " _ 88W 

du \du ) \du" ) 88u 


(7a) 

(7b) 


along with static boundary conditions that are naturally satisfied due to the 


condition on the displacements v and u. All derivative terms in Eq. (7b) 
reported in Appendix [A] 


27T periodicity 
are explicitly 


B. Asymptotic expansion and reactive force of the substrate 

We seek a solution of Eq. (J7|) as an expansion of the form 

v — Vo + eA sin ( m6 ) + O (e 2 ) , (8a) 

u = eB cos ( md ) + O (e 2 ) , (8b) 

where (uo ; 0) corresponds to a radial pre-buckling deformation and e(Asin (md ), B cos (md)) 
represents an instability of azimuthal wavenumber m. From the requirement of ‘In periodic 
functions v and u, m has to be an integer. Also, we consider m > 1 since m = 1 corresponds 
to a solid body translation of the ring. The instability displacement held has amplitudes eA 
and eB, where £ is a small parameter. 

Before substituting Eq. Q into Eq. ([7]) and solving at each order in e, we hrst need to 
determine the work, 8W , done by the reaction force, F = — ((rej_ + rey), that the substrate 
exerts on the ring. The normal stress a and the tangential stress r at the interface between 
the ring and the substrate are computed by solving the corresponding two-dimensional 
elasticity problem using an Airy function, with pressure P at r = R$ and the displacement 
held given by Eq. ([8]), at r = R. 
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It is worth to note that in our computation, we enforce the continuity of the displacement 
and stress fields at the interface ring/substrate, similarly to |H)]. This is a main difference 
with previous studies of wrinkling surfaces, e.g. in [33], who assumed zero shear stress and 
continuity of the normal stress. Such an approach yields a stretching energy in the ring 
much larger than the bending energy, in contradiction to what is expected in the wrinkling 
of a thin him m Our explicit solution for the boundary value problem of the substrate is 
reported in [B] 

In short, we find that the stresses a and r at the interface are 


= = k 0 i>o + keA sin (md) — yp, 
E F 

T 

= = peB sin (mO), 

E F 

with the following governing parameters 


k 0 = 
k = 


(1 - i/g) (1 - l 2 ) 1 

1-2 v s +F r 


(1 - o?) 


Sa + S>b 


1 

r 


k = 

7 = 


(1 - vs) 
2 



W 2 (1 - Us) 
l-2v s + P 2 ' 


1 

r 


(9a) 

(9b) 

( 10 a) 

( 10 b) 

( 10 c) 

(lOd) 


Here, Ay, k, and // are the dimensionless pre-wrinkling, wrinkling and shear stiffnesses of 
the substrate, respectively. The coefficient 7 quantifies the effect of pressure, p, on the ring 


through the term —yp in Eq. (9a). The value of y is smaller than one, reflecting the fact that 
the transmitted pressure decreases within the substrate. There are however two exceptions 
for which y = 1 : the limit where there is no substrate, /3 —» 1 , and the limit of a perfectly 
incompressible substrate, v$ —> 0.5, in which the volumetric pressure remains constant 


throughout the substrate. The quantities Sa, Sp and Tp used in Eq. ( 10 ) are functions 
of both the cavity size /3 and the azimuthal wavenumber m, and their full expressions are 
reported in |Bj 

From the solution of the linear elasticity problem for the substrate (see 0 , the reaction 
force, F. has a constant direction along e r . However, based on physical intuition, one would 
expect the restoring force to change its direction as the ring deforms. Accounting for this 
scenario in full would have required solving the elasticity problem for the substrate, with 
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nonlinear kinematics. Here, for simplicity, we assume that the magnitude of the restoring 
force is given by the solution of the linear elasticity problem, while the direction of the force 
is given by the vectors e_L and ey. Depending on which configuration of the ring is used 
(initial or deformed), two cases need to be considered to define the normal and tangential 
vectors e_i_ and ey, respectively. If we define these vectors from the initial configuration, then 
(e_i_,ey) = (e r ,ee), so that the reaction force F keeps a constant direction while the ring 
deforms. If we use the deformed configuration, then F is modeled as a follower force whose 
direction changes with the ring deformation. In this case, (e_|_, ey) = (N. T), where N is the 
inward normal vector orthogonal to the ring center-line, in its deformed state. In order to 
take both of these options into account in the same model, we define a parameter, y, which 
can take the values y = 1 or y = 0, when either the undeformed or deformed configurations 
are used. Thus, we express the reaction force and its elementary work per unit length of the 
ring as 

F = - [(cre r + re e ) y + (ctN + rT) (y - 1)], (11a) 

5W — FR (Sve r + 5ue e ). (lib) 

where • is the Euclidean dot product. 


1. The principal solution: order 0 in e 


We now proceed to obtain the principal solution and the critical instability modes by 
substituing Eq. (J8| into Eq. ((tJ) and solving the resulting equations for each order of e. At 
order 0, Eq. (J 3 J for the stretching deformation 7] yields r/o = Vq = cr 0 /Ep and the linear 


approximation of the Euler-Lagrange Eq. (7a) 

IP _ 

h + h 3 /12 + ko Ep 


Vo = 


( 12 ) 


relates the dimensionless pressure, p, and the dimensionless hoop stress, ao/Ep, in the ring. 
Note that Eq. © is automatically satisfied at order 0 in e. In the absence of a substrate, 
the ratio a 0 /P reads 


Hhi&) = l +0(h) ’ 

and in the limit of the ring that is much stiffer than the substrate we obtain 


lim 

£—>■00 




(13) 


( 14 ) 
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so that < lim (^) < lim (^). For vg — 0.5, we have 7 = 1 through the dehnition 

r £—^00 v ^ ' / 3 -U v ^ ' 


m 


Eq. (lOd) and both limits in Eqs. (d3h and (14) are equal. Thus, the hoop stress in a very 


stiff ring lying on a soft incompressible substrate is the same as the hoop stress in a ring 
with no substrate, which serves as a verification of the rationale thus far. 

2. The instability: order 1 in e 


At order 1 in e, the Euler-Lagrange Eqs. ([7]) write 

(dip + ai) A + [hip + b\ j B = 0, (15a) 

(< a 2 p + ajj) A + (b 2 p + b^j B = 0, (15b) 

where a*, a), bi and bi are functions of the azimuthal wavenumber m, the stiffness ratio £, the 
cavity size (3 and the substrate Poisson’s ratio ug, reported in full in[Cj The linear system 


in Eq. (15) has a nontrivial solution when its determinant vanishes, which for a given value 


of m occurs at the dimensionless pressure 

P m —dib 2 + bid 2 

Pm — == — ~ ~ _ ~ - 

Bf dib-2 + dib-2 - bid 2 - b\d 2 

The critical azimuthal wavenumber and the dimensionless critical pressure can now be ob- 


( 16 ) 


tallied from Eq. (16) by minimizing over all possible values of m 


[m c ,p c }= min (p m ) . 

m= 2,3,... 


( 17 ) 


For the general case, this minimization must be performed numerically. However, in the 


absence of the substrate, we obtain that m c = 2 and Eq. (16) simplifies to 

1 


lim (p c ) = lim (p 2 ) - - 
/3—>•! /3-U 4 — X 


- 


;h 5 + O ( h ? ) 


(18) 


(4 — xY 

The term of order h 3 corresponds to the classical dimensionless critical pressure for a ring 
with no substrate PSISM]. The term in IP is a correction also reported by [50] • In the 
limit of a ring much stiffer than the substrate, we also have m c = 2 , with 

1 (-1 + 2 v s -p 2 )h 3 2 (—1 + 2u s — (3 2 ) h 5 


lim ( p c ) = lim (p 2 ) = — , „ 

C^oo 5^00 2 [3 2 (-1 + u s ) (4 - x) (3 2 (-1 + ug) (4 - xY 


+ O ( h 7 ) , (19) 


so that lim ( p c ) < lim ( p c ) < 0. Again, for ug = 0.5, we observe that the limits equal, 

£->■00 / 3 -U 

corroborating the physical intuition that as the stiffness ratio £ increases, the effect of the 
substrate becomes negligible. 
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C. Critical stress, wavelength and comparison to wrinkling on a planar substrate 


To determine the effect of curvature on the instability, we compare the critical hoop stress, 


cr c , in the ring, given by Eq. (12) with p = p c , as well as the wavelength A c = 27 t R/m c of the 
wrinkling mode, against their counterparts for an initially planar him on a plane substrate 
of infinite thickness m 


/ ^*\ 2/3 — 


E 


E s 


OPlane . I 3 —— 1 


Er 


Ef ( 4(1 — vs)~ 3 
3-4 u s e 


2/3 


'Plane 


= 2irH - 


n e f x 1/3 


V3 E s 


= 2*ml 3 ~ 4 ^ 


3 4(1 - vsY 


1/3 


(20a) 

(20b) 


where Es* = 4As(l — vs ) 2 /{3 — 4z/g) is the effective stiffness of the substrate. Together, 


Eqs. (|20j) and (jl2j) yield 
4 7Pc 


Or 


OPlane ^ T j 12 T /lq 

If 3 - Au s 


^Plane hm c V 12(1 - v s ) 



- 2/3 


(21a) 

(21b) 


and, in the limit of absence of the substrate, Eq. (]21a|) reduces to 

Or. \ i .. , . i ,. , , 4 / 12(1 - v s ) 


lim 

\<7pi 


- 2/3 


OPlane 


lim (of) = 


-lim (oo) = . . , . ^ 

Opiane Z 3 ^ 1 4-^1 (3 - 4l/g) f 


h 2 + 0(/r 4 ) , 
( 22 ) 


where i s the Poisson’s ratio of the substrate in the planar case. 


In § VD we shall make use of Eqs. (21) and Eq. (22) to further quantify the sensitivity 
of the stress and wavelength to the curvature of the substrate. 


IV. NUMERICAL SIMULATIONS 

In § |V| we will contrast the predictions from the above analysis with the results of a 
series of finite element simulations performed using the commercial package Abaqus, with 
the BUCKLE analysis, which provides the buckling load and the corresponding eigenmodes. 

As in the analytical study, the cylindrical structure was modeled as an annulus of a 
soft substrate with a stiffer thin him adhered to its exterior, under plane strain conditions. 
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The bonding between the substrate and the film is assumed to be perfect, such that both 
share nodes. Negative pressure is applied on the interior surface, to model the effect of the 
pressure differential between the inner surface of the annulus (the cavity) and the exterior 
of the system. The buckling analysis provides the value of the critical pressure, p c , as well 
as the corresponding critical mode. Rigid body motions are removed by constraining the 
displacement of two points on the film. 

Both substrate and film are modeled as incompressible linearly elastic materials, = 
vp = 0.5. These results were compared to additional simulations using a Neo-Hookean 
model but no difference was observed given the low values of strain involved. The substrate 
was modeled using quadrilateral plane strain elements. Due to the incompressibility, the 
corresponding hybrid element, CPE4H, was used. The film was modeled using B21 beam 
elements. In order to account for the effect of plane strain, the stiffness of the beam is 
defined as E F / (1 — u F ). 

All of the results presented were obtained using 1000 elements in the circumferential 
direction, and 150 (R — Rq) elements in the radial direction. The mesh size was validated 
with a convergence analysis. By way of example, differences less than 0.5% in critical pressure 
were obtained when comparing results for a mesh with twice the elements in each direction, 
even in the cases of wrinkling with the shorter wavelengths. The deviations between the 
two meshes were, however, larger (~ 5%) for the case when the critical buckling mode is a 
Biot instability, due to the infinite number of wavelengths associated to the same buckling 
mode. However, as stated in § [TTJ this mode shall not be studied in detail as we focus on 
the wrinkling and global buckling loads. 

To test the validity of using beam theory to describe the ring in our problem, we also 
performed numerical simulations with 2D solid elements (CPE4H, the same used for the 
substrate) and found excellent agreement with the simulations using beam elements. How¬ 
ever, in order to achieve such agreement, the mesh needs to be greatly refined, resulting in 
a significant increase in computational cost. Attempts to use a mesh size similar to that of 
our previous simulations showed clear disagreement. Given the excellent agreement between 
the two versions - either using CPE4H (with a fine mesh) or B21 for the ring - as well as 
the significantly lower computational cost of the B21 elements, we have decided to use beam 
elements for the ring in our analysis. 
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V. RESULTS 


Having introduced our analytical and numerical methods, we now present the results of 
a systematic exploration of the mechanical response of our system for different geometric 
and material parameters. Throughout, we provide a direct comparison between analytical 
results and numerical simulations, finding good agreement. A few instances of discrepancy 
will also be discussed. 

For the geometric parameters, we have varied the cavity size, f3 = Ro/R, and the dimen¬ 
sionless thickness, h = H/R. Three representative values were chosen for /3 (depicted in the 
insets of Fig. [3^-c): a small cavity, /3 = 0.2; a cavity with size half of the external radius, 
/3 = 0.5; and a large cavity, ft = 0.8. Moreover, h was varied in the range 10 -3 to 10 _1 . 
This parameter has two different physical interpretations. On one hand, for a substrate with 
given curvature, i.e. fixed R , increasing h is equivalent to increasing the thickness of the 
ring. On the other hand, for a ring of given thickness H , the value of h decreases with the 
curvature, 1 / R. 

For the material properties, we have considered values for the stiffness ratio between the 
film and the substrate, £ = Ep/Eg, spanning over five orders of magnitude, from 10 2 to 10'. 
The substrate is taken to be incompressible, vg — 0-5, since most of the relevant experiments 
that have motivated our study mmmm use nearly incompressible elastomeric sub¬ 
strates. Moreover, it is important to note that, even if the analytical model has been derived 
assuming a general value of ug, it is expected to be less accurate for increasing deviations 
from incompressibility. This was analyzed by [30], who showed that for the wrinkling of 
plates deviations when v s = 0.3 are of just a few percent. 

In the presentation of our results, we first consider the effect of h and £ on the hoop-stress 
and on the critical pressure. Then, we rationalize the transition from wrinkling to global 
buckling, shown in Fig. [l] We finally compare the critical stress and wavelength of the 
wrinkling mode to their planar substrate counterparts and discuss the effect of curvature. 


A. Hoop stress prior to wrinkling 

In Fig. [3j we plot the hoop stress in the ring, normalized by the pressure, (Jo/P, as a 
function of the dimensionless thickness, h = H/R. We find that the hoop stress decreases 
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monotonically with the dimensionless thickness, and increases with both the stiffness ratio 
and the cavity size. When the ring is much stiffer than the substrate (i.e. £ —y oo), the hoop 


stress scales as cr 0 /P ~ ( H/R ) _1 , with a prefactor given by Eq. (14). For a given thickness 
of the ring and a given cavity size, the hoop stress in the ring decreases as the curvature of 
the substrate increases. This observation is consistent with the classic result for the hoop 
stress, do = PR/H, for a depressurized thin-walled cylindrical pressure vessel. This result 


can be recovered from Eq. (12) by taking 7 = 1, ko = 0 and performing a Taylor expansion 
in h, about 0 . 


B. Critical pressure 


In Fig. |4| we plot the dimensionless critical pressure, \p c \ = \P C \/Ep, as a function of 
the dimensionless thickness, h, and observe two different regimes. For low values of h, the 
ring wrinkles with an azimuthal wavenumber m c 2 , with a critical pressure \p c \ that 
increases with h. When h reaches a threshold value, h*, the wavenumber decreases suddenly 
to m c = 2 ; the nature of the instability changes from a wrinkling to a global buckling mode. 


In this regime, there is a clear asymptote as £ —> 00 , given by Eq. (19), which corresponds 
to the buckling of a ring with no substrate. 

Depending on the model used for the reaction force of the substrate (constant direction 
force, x = 1, or follower force, x — 0, see § III B), Eq. yields two different analytical 
predictions, shown in Fig. [4] as a solid line for x = 0, and a dashed line for x = 1- We 
note that the agreement between FEM simulations and analytical predictions is superior for 
X = 0, in particular for the global mode. In other words, it is better to model the reaction 
force of the substrate as a pressure held which remains normal to the ring center-line C as it 
deforms, than as a force with constant direction. Thus, from now on, all analytical results 
will be only presented for x = 0 . 

Despite the overall good agreement between the analytical and numerical results, there 
are noticeable discrepancies in the global modes for low values of the stiffness ratio. The 
reason is that, for £ = 10 2 , the strains in the substrate can become significant, such that the 
assumptions for our linear theory are no longer valid (see 0 - A particularly extreme case 
is the appearance of Biot modes for /3 = 0.2 and £ = 10 2 . These modes, although possible, 
only occur in a small region of our parameter space and, as mentioned above, are beyond 
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FIG. 3. Pre-instability hoop stress, <ro, in the ring, normalized by the pressure, P, as a function 
of the dimensionless ring thickness, h = H/R. Cavity sizes are: (a) (3 = 0.2, (b) (3 = 0.5 and 
(c) /3 = 0.8. The Poisson’s ratio of the substrate is = 0.5. Analytical predictions are given by 


Eq. (12) as solid lines and FEM results are plotted as data points. The legend (bottom right) is 
common to all three plots. 


the scope of this work. 

In short, our model exhibits limitations if the ring and the substrate have comparable 
stiffness, or when there is a Biot mode (small cavity). Apart from this extreme combi¬ 
nation of parameters, rarely observed in experimental configurations, the model performs 
successfully. 
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FIG. 4. Dimensionless critical pressure \p c \ = \P C \/Ep versus the dimensionless ring thickness 
h = H/R, for several stiffness ratios £ = Ep/Eg and cavity sizes (3: (a) (3 = 0.2, (b) (3 = 0.5 and 


(c) (3 = 0.8. Theory (lines) is given by Eq. (17). The Poisson’s ratio of the substrate is 1/5 = 0.5. 
The legend (bottom right) is common to all three plots. 


C. Phase diagram 

We proceed by focusing on the transition from wrinkling to global buckling, towards 
first constructing a phase diagram in the (£, h) parameter space and then quantifying the 
dependence of the boundary, h*(£), between the two modes on the size of the cavity, (3. 
Given that we do not have a closed form expression for h*, we use a numerical method 
which tracks any jump from m c = 2 to m* c > 2 in Eq. 0 , when h is decreased. 

In Fig. |5](a), we plot h* as a function of £, for /3 = 0.5. The boundary between modes 
is consistent with a power-law, h* ~ £~ 3 , which divides the phase diagram into a wrinkling 
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FIG. 5. (a) Phase diagram in the (£, h) parameter space, showing the transition boundary, h*, 
from wrinkling to global buckling, for (3 = 0.5. (b) Phase diagram for j3 = 0.2 to (3 = 0.8. Solid 
lines are analytical predictions and dotted lines correspond to the numerical fit £ = f 2 h 21 . (c) 
Fitting coefficients, f \ and / 2 , as functions of (3. (d) Azimuthal wavenumber of the wrinkling mode 
for h = h*, as a function of /3. The Poisson’s ratio of the substrate is v$ = 0.5. 

domain (h < h*) and a global buckling domain (h > h *, shaded region). Again, there is 
good agreement between numerical and analytical solutions. 

In Fig. |(b), we extend the phase diagram to cavity sizes from /3 = 0.2 to [3 — 0.8. We 
predict that the power-law with exponent —3, mentioned above for /3 = 0.5, is still valid as 
(3 increases, even if there are some deviations towards the higher values. To quantify the 
appropriateness of the —3 power-law, we fit a curve of the form £ = / 2 /P 1 to the analytically 
calculated boundaries. In Fig. Kc) we plot the fitting parameters f± and /2 as a function of 
f3. We find that the exponent is fi ~ —3 for (3 < 0.7, and then decreases rapidly to a value 
of fi ~ —4.25 for (3 = 0.8. The decrease of the prefactor f 2 is more pronounced and reflects 
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the fact that h* decreases as /? increases. 

The evolution of m* (i.e. azimuthal wavenumber of the wrinkling mode for h = h*) as a 
function of (3 and £ is presented in Fig. |5](d). We observe that m* increases with f3, while 
it is nearly insensitive to £. For example, for /3 = 0.5, we find that m* = 9 for £ > 10 2 and 
7Ti* = 10 for £ = 10 2 . Predictions for £ > 10 2 are in relatively good agreement with finite 
element simulations, which have shown that m* = 10. For £ = 10 2 , FEM simulations yield 
77i* = 15, pointing out, once again, the limitation of our analytical approach in the limit 
when the stiffness of the ring and the substrate become comparable. 


D. Comparison of wrinkling in our curved system with that on infinite planar 
substrate 


Thus far, we have shown that the ring may wrinkle or buckle globally, depending on 
the curvature and stiffness of the substrate. We now focus on the wrinkling mode of the 
ring, with the aim of comparing the critical stress and wavelength of our curved system 
to their counterparts for a infinite planar substrate. We shall center our discussion of this 


comparison for /3 = 0.5. In D we report the results for /3 = 0.2 and /3 = 0.8, which are 
qualitatively similar. 


In Fig. [ 6 |a), we plot cr c /api ane given by Eq. (21a), as a function of h. For h < h* 
(wrinkling domain), we find that <J c /crpi ane ~ 1, in agreement with the FEM simulations. 
To further quantify how close to unity is this ratio, in Fig. | 6 ](b) we plot the dimensionless 
critical normal force a c h/Ep as a function of the planar result, opi ane h/Ep. We obtain a 
line with unit slope, indicating that the substrate curvature has no significative effect on the 


critical stress for wrinkling. From this observation, and using Eq. (12) with (Tq = crpi ane , we 


write the following approximation for the dimensionless critical pressure, 
h T />'() (T plane. 


\Pc I = 


+ O (h 3 ) , 


(23) 


7 E F 

with k 0 and 7 given by Eq. (10). Finally, back to Fig. [ 6 ](a) , for h > h* (global buckling 
domain), we find that cr c /api arLe first decreases with h, then reaches a local minimum and 
eventually increases as a power-law with slope 2 , towards the asymptotic limit given by 


Eq. (22). The evolution of <J c /<Jpi an e in the global buckling domain is also reproduced well 
by FEM simulations. 
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FIG. 6. (a) Critical stress a c normalized by its planar substrate counterpart, crpiane > as a function 
of the dimensionless ring thickness h = H/R. (b) Dimensionless critical normal force versus planar 


substrate counterpart. Analytical predictions (solid lines) are given by Eq. (21a) and FEM results 
are shown as data points. Insets are sketches of the wrinkling and global buckling modes. The 
cavity size is (5 = 0.5 and the Poisson’s ratio of the substrate is vs = 0.5. The legend is common 
to both plots. 


We now investigate the effect of curvature on the wavelength of the instability mode. 
In Fig. [7|a), we plot X c /Xpi ane , given by Eq. (21b), as a function of h. Focusing on the 
wrinkling domain, we find that A c /A pi an e ~ 1, in agreement with FEM simulations. In 
Fig. 0b), once again, we quantify how close to unity this ratio is by plotting the critical 
wavelength A c as a function of the planar result, Xpiane■ To Erst approximation, noting the 
large dynamic range (at least two orders of magnitude) in both axes of the plot, we find a 
line with unit slope that passes through the origin, suggesting that the substrate curvature 
has no significative effect on the wavelength of the wrinkling mode. The only deviations 
arise from the discrete nature of the wavenumber, since the geometry of the ring enforces 2n 
periodic wrinkling modes. To highlight this phenomenon, in Fig. [7](c) , we plot the critical 
azimuthal wavenumber m c , given by Eq. 0. as a function of h and we superimpose the 
planar substrate result (dashed line), mpi an e = 2irR/Xpi ane - We find that m c is a decreasing 
stair function of h. The deviation in wavelength between the curved and planar cases scales 
as (|A C — Xpi ane \)/Xpiane ~ 1 /m, which is maximum for m*. As shown previously in Fig. 
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FIG. 7. (a) Critical wavelength, A c , normalized by its planar substrate counterpart, Ap/ ane , as 
a function of the dimensionless ring thickness h = H/R. Analytical predictions (solid lines) are 


given by Eq. (21b). (b) Zoom of the wrinkling domain by plotting A c versus A plane- ( c ) Azimuthal 
wavenumber m c of the wrinkling mode as a function of h. Analytical predictions (solid lines) are 
given by Eq. ©• Dotted lines show m c for the planar substrate case. Insets on (a) and (b) 
are sketches of the wrinkling and global buckling modes. Inset on (c) is a zoom in showing the 
difference in m c for the curved and planar substrates. The cavity size is (3 = 0.5 and the Poisson’s 
ratio of the substrate is us = 0.5. FEM results shown as data points. The legend (bottom right) 
is common to all three plots. 


id). m* increases with /5, hence the deviation (|A C — Xpiane |) /^piane is maximum for a small 
cavity size. By way of example, for f3 = 0.5, Fig. |5](d) indicates that m* = 9, leading to a 
maximum deviation (|A C — Xpi ane \) / Xpi ane ~ 10%, that rapidly decreases as m increases. 
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VI. CONCLUSION 


We have considered the two-dimensional problem of a ring bound to an elastic substrate 
which contains a cavity that is depressurized. An energy formulation was used to derive 
the Eulcr-Lagrange equations that govern the equilibrium of the ring, and solved them via 
an asymptotic expansion. As an improvement to previous results in the literature, our 
analytical approach accounts for the effect of curvature in modeling the reaction force of the 
substrate. These analytical results were compared with numerical simulations. 

We Erst studied the principal solution, obtaining an expression for the hoop stress in 
the ring as a function of the applied pressure. We then performed a stability analysis of 
the problem to determine the critical pressure, P c , and the corresponding instability mode. 
Depending on the dimensionless thickness and stiffness ratio (h and £) we have identified 
two different regimes: local wrinkling of the ring, and global buckling of the structure. The 
boundary between both regions of instability was described via a detailed phase diagram, 
which quantifies the value of h and £ at which the transition between instabilities occurs 
and takes into account the cavity size, f3. Our results can be used as a design guideline to 
target a desired mode. Finally, we have shown that the critical stress for wrinkling and the 
resulting wavelength do not depend significantly on the curvature of the substrate. However, 
curvature imposes a discretization of the wrinkling wavelength due to the periodic closing 
conditions of the ring. 

Our study focused on a 2D curved system which exhibits instability modes analogous to 
the cylindrical pattern found for uniaxial compression of a him on a hat infinite substrate. 
Considering more complex loading conditions ( e.g. also introducing axial loading) or shells 
with non-zero Gaussian curvature, should lead to more complex patterns that deserve to be 
investigated further. Having validated the FEM analysis, as well as carefully considering 
the elastic response of the curved substrate, extending the study to these other scenarios 
and addressing the issue of pattern selection is an exciting avenue for future research. 
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Appendix A: Derivative terms in Euler-Lagrange equation 


The terms in the Euler-Lagrange equation for the equilibrium of the film, Eq. (7b), are 


as 

dv 


du 






~- (-2m 2 + (2 - 4m") v + 2m" - v' 2 + 2mm" - 4u'v" - 2 v'u" + 2 u' 2 + u 2 
^ (m 2 + (2 + 2 u') v + u' 2 + 2 u' + (m - m') 2 ) (1 + v + u') (Ala) 

^ (-2m 2 + (2 - 4v") v + 2v" - v' 2 + 2mm" - 4 u'v" - 2 v'u" + 2 u' 2 + m 2 ) (u + u") 

^ (m 2 — 2uv' + v' 2 + (v + 2 + m') (m + m') j (m — v') (Alb) 

_ I2 (( -2v + 1 ~ 3V " ~ U ‘ W ) ^ + ^ _ 2v '" + “)“'- 3 m"m" + (! - 2m) m'" + «/) (m' + 

^ ((—4m + 2 - 4u') v" - 2v 2 + 2v - 2v'u" - v' 2 + 2mm" + 2 u' 2 + m 2 ) (m" + u '") 

(u - v') ((v + 1 + v") v' + (u + u") u' + ( 1 + v) u" - uv") 

^ {v! — v") (m' 2 + (2m + 2) u! + v 2 + 2v + (u — m') 2 j (Ale) 

— (m' — v") ((—2m + 1 — 3v" — u'") v' + (3 u" — 2 v'" + u) u' — 3 u"v" + (1 — 2v) v"' + uu 1, 
6 

^ ((2 u - 2v') u" + (-4m + 2-4 u') v" - 2v 2 + 2v + 2 u' 2 - v' 2 + m 2 ) (-m'" + m") 



((« + m") m' + (m' + m") m + m" + (1 + v") v' - uv") (1 + V + u') 

^ (m' + u") (m ' 2 + (2v + 2) u' + v 2 + 2v + (u — m') 2 j (Aid) 

~ y { v ~ 1 +u ') (~ 3v " 2 + (i _ 2v - Au " r) v " + u ' 2+ ^~ 2v "" +Au " r) u>+{i - v "") 

J (v - i + m') (3m" 2 + (—5m'" + m) m" - 2m' 2 + (-5m'" - m"") v + ««'j 

((—2u + 1 - 3v" - u m ) v' + (-2 u' + l-2v) v m + (3m" + u) u' + uu'" - 3u"v") (v' + u 
o 

^ ((-4u + 2 - 4 u') v" - 2v 2 + 2v - 2v'u" - v' 2 + 2 uu" + 2 u' 2 + u 2 ) (w" + u'") 

(Ale) 


2- (m - v') (-3v" 2 + (1 - 2v - 4 n'") w" - 2^ /2 + (-5t/" - m / w ) m' + (-2m' + 1 - 2u) m""' 
lz V / 

^ (m - m') ((«"" + m") m - 5m"m'" + 3m" 2 + 4u!u!" + m' 2 ) 

— (m' - v") ((-2m + 1 - 3m" - u'") v' + (3m" - 2m"' + u) u' - 3m"m" + (1 - 2m) m '" + mm'" 
6 

^ ((2m - 2m') m" + (-4m ?- 5 2 - 4m') m " - 2m 2 + 2m + 2m' 2 - v' 2 + m 2 ) (- m '" + u") 


( Alf'i 





and 


ddW 

dSv 


ddW 

ddu 


1 2 r (-1 + x) (-1 + u' + v) (u - v')- a ((-1 + x) v' 2 - 2 u (-1 + x) v' + 2 + (-1 + x) u 2 ) 

2 h 


(A2a) 

1 -2 a (-1 + v! + v) (u - v') (-1 + x) ~ T ((-1 + X) v' 2 - 2 u (-1 + x) v 1 + 2 + (-1 + x) u 2 ) 

2 h 


(A2b) 


Appendix B: Response of the substrate 


In this appendix we consider the boundary value problem of the substrate, subjected 
to the ring displacement R (v (0) ,u(0)), v and u given by Eq. (J8]) , at the interface r = R 
between the ring and the substrate and to the pressure P at r = R 0 = f3R. The substrate 
is assumed to be in a state of plane strain. We note a = 3 — 4and introduce the shear 
modulus G = Es (1 — vs)/ 2. The components of the stress in the substrate are represented 
by oy r , cr r o, and the displacement field is ( U r , Ug). The 2D problem of elasticity is solved by 
finding an Airy function of the form |52] , 

(j) m (r, 9) = B ir 2 + B 2 In (r) + £ ( A x r m+2 + A 2 r~ m+2 + A 3 r m + A A r~ m ) sin (mO) , (Bl) 
where A, and Bi are unknown constants determined by the boundary conditions 


a rr (i? 0 , 6) = —P, (B2a) 

a r g (Ro, 6) = 0, (B2b) 

U r (P, 9) = Rv ( 9 ), (B2c) 

U e (P, 9) = Ru ( 9 ), (B2d) 


which are assumed to apply at r = P 0 and r = P. The first two equations stand for the 
continuity of the stress at the boundary of the cavity, whereas the last two stand for the 
continuity of the displacement at the interface between the substrate and the ring. 


The stress and displacement fields resulting from the Airy function Eq. (Bl) are [53 


a rr (r, 9) = 2Bi + 



—A 1 (m + 1) (m — 2) r m — A 2 (m + 2) (m — 1) r m 
— A 3 m (m — 1) r m ~ 2 — A 4 m (m + 1) r -m ~ 2 


sin ( m9 ) , 
(B3a) 
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cr re {r, 6) = 


—Aim (m + 1) r m + A 2 m (m — 1) r m 
—A 3 m {m — 1) r m ~ 2 + A^rn (m + 1) r ~ m ~ 2 


cos (m6), 


(B3b) 


= ia 


Bo I Ai(a — m — 1 ) r m+1 + A 2 (a + m — 1 ) r -m+1 

-Bi (a - 1) r -- + 

r ' —A 3 mr m ~ 1 + 


sin (■ m6 ) 
(B3c) 


1 / -Ai(a + m + l)r m+1 + A 2 (a-m + l)r m+1 

1 |cos(m0). (B3d) 

2G l —A 3 mr m 1 — A^mr m 1 


Applying the boundary conditions Eq. (B2) to Eq. (B3) yields a linear system for Ai and 
Bi, with solution 

1 


Bi = 


2 (/3 2 + 1 - 2 u s ) 


[E s v o (1 - u s ) - P/3 2 ] , 


(B4a) 


Bo = 


R P fTW 


/3 2 + 1 - 2i/ 5 


[B 5 u 0 (1 - u s ) + B (1 - 2i/ s )] , 


(B4b) 


A 1 = 


(m — 1) /3 2m+2 + (a — m + 1) /3 2m + /3 2 


2 (— m 2 + 1) f3 2m+2 + m 2 /3 2m + 4 + /3 4m+2 a + {m 2 — 1 + a 2 ) /3 2m + ct/3 2 
(m - 1) / 3 2m+2 + (—m - a + 1) ,5 2m - /3 2 


2 (—m 2 + 1) /3 2m+2 + m 2 P 2m+A + /3 4m+2 a + {m 2 — 1 + a 2 ) /3 2m + a/3 2 


GR~ m A 


GR~ m B , 


(B4c) 


Ao — 


(m + 1) /3 2m+2 + /3 4m+2 + /3 2m (a + m + 1) 


2 (—m 2 + 1) /3 2m + 2 + m 2 /3 2m+4 + /3 4m+2 a + (m 2 — 1 + a 2 ) /3 2m + a/3 2 
(m + 1) /3 2m+2 + /3 4m+2 + /3 2m (a - m - 1) 


2 {—m 2 + 1) /3 2m+2 + m 2 P 2m+A + /3 4m+2 a + (m 2 — 1 + a 2 ) /3 2m + a/3 : 


GR m A 


-GR m B , 


(B4d) 


^ 2 ((a - m + 1) (m + 1) ,d 2m + m 2 /3' 2m+2 + a + 1 + m) CRr m+2 4 

(2 (—m 2 + 1) /3 2m + 2 + m 2 /3 2m+A + /3 4m+2 a + ( m 2 — 1 + a 2 ) /3 2m + a/3 2 ) m 

/3 2 {m 2 p 2m+2 - 1 + (—m - a + 1) (m + 1) /3 2m - m + a) „ m+2 

(2 (—m 2 + 1) /3 2m + 2 + m 2 /3 2m+4 + /3 4m+2 a + (m 2 — 1 + a 2 ) /3 2m + a/3 2 ) m 
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(B4e) 


(- m 2 j3 2m+2 + (m - 1) (a + m + 1) / 3 2m - (a - m + 1 ) /3 4m ) /3 2 m+2 . 

(2 (—m 2 + 1) /3 2m + 2 + m 2 /3 2m+4 + /3 4m+2 a + (m 2 — 1 + a 2 ) /3 2m + a/3 2 ) m 

(m 2 /3 2?n+2 + (m — 1) (a — m — 1) /3 2m — (—m — a + 1) /3 4m ) /3 2 GR m+2 B 
(2 (—m 2 + 1) /3 2m + 2 + m 2 j3 2m+i + /3 4m+2 a + (m 2 — 1 + a 2 ) /3 2m + a/3 2 ) m 

(B4f) 


Substituting A, and B, into Eq. (B3) yields the stress at the interface r = R 


a rr (R, 9) = A'o-R^o + KReA sin ( mO ) — 7 P, 

cr r g ( R, 6) = MReB cos ( mO ), 


where 


and 


K _ 1 (1 — u s) (1 ~ /^ 2 ) 

0 Es R l-2v s + (3 2 ' 


K = E 


1 (i - v s ) 


' R 


M = £ s - 


-w 1 (! - ”s) 


7 = 


R 2 
2/3 2 (1 - u s ) 
1-2 vs + P 2 ' 


Sa + <Sb — 

S b ^+T b 


i/jSa = 2 (— m 2 (/3 4 + 2 + a) + (3 + a) (m 2 — l) /3 2 + 2 (a + 1)) /3 
— (((a + 1) m + a — 1) /3 4m — (a + 1) m + a — l) /3 2 , 


2m 


(B5a) 

(B5b) 

(B 6 a) 

(B 6 b) 

(B 6 c) 

(B 6 d) 

(B7a) 


ipS B = —2 m (— m 2 (/3 4 + l) + 2 (m 2 — l) (3 2 + 1 + a) /3 2m 
+ (((a — 1) m + a + 1) /3 4m + (a — 1) m — a — l) /3 2 , 


(B7b) 


^T b = 2 (—m 2 /3 4 + (1 - a) (m 2 - l) + am 2 ) /3 2m 

— (((a + 1) m + a — 1) /3 4m — (a + 1) m + a — l) /3 2 , 


(B7c) 


= (m 2 f3 A - (m 2 - 1) (2 /3 2 - 1) + a 2 ) /3 2m + a/3 2 (l + /3 4m ) . 


(B7d) 









The dimensionless stress at the interface and the dimensionless stiffness parameters are 


obtained from Eq. (B5), by dividing it with E F , 


a 


= = k 0 v 0 + keA sin(m0) — qp, 
E F 


T 


= = iieB sin(m0), 
E F 


K 0 R _ ( 1 - u s ) (1 - f3 2 ) 1 
1-2 u s + (3 2 £’ 


fcn= = = 




aj t 

MR (l-u s )f A ^ 1 

H = -= = --- [S B +T B -, 


Er 


B 




(B8a) 

(B8b) 

(B8c) 

(B8d) 

(B8e) 


as indicated in Eq. (10) 


We note that in the case of an inextensible wrinkling mode, A/B = m, the stiffness K 
simplifies to K given by 


K = E,s4 (1 0 US) ( S a +S b - 
R 2 V m 


= E S 1 


mR 


2(1 - u s ) ( f3~ 2m - f3 2m + 2m (1 - (3~ 2 )) (m 2 - 1) 
tf- 1 - j3) 2 (m 2 - 1) + (/3— - (3 - 4 u s ) + + p) 


2 i 


(B9) 


which, for a substrate with no cavity, leads to the limiting case 


Mm (K) = 

/?—»o V J mR 3 — 4 us 


(BIO) 


and, for an infinite plane substrate, yields 


lim 

R —^oo 


( lim 

\/ 3^0 



B Plane 


— 4(1 - U S ) 2 7T 
S 3-4 u s A’ 


in agreement with [27] . 
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Appendix C: Terms of the linear stability analysis 


The terms a,- 


Zj, bi and b, that appear in Eq. (16) are 


a i 


a i 


bi 

h 


Oj 2 


b 2 

b 2 


—A^h f f m 4 + - — — m 2 ^ h 2 — 9 — 3m 2 ) -- -J- -—, 

^ VV 2 4/ /z + /z 3 /12+fc 0 ’ 

h ((m 2 - l) 2 h 2 + 12 ) £ + 6 (1 - u s ) S A , 

—A^h ( -h 2 m 3 + 12 m ) -- —p -—, 

S V 4 ) h + h?/l2 + ko 


-12 £hm + 6 (1 - u s ) S B , 

- ( h 3 m 2 + A8I1 + 12 k 0 (1 - *)) m£ + ^ + ^ + 12y (1 - X ) m£, 
—12hm£ + 6 (1 - u s ) S B , 

(36 hm 2 + h 3 + 12 h + 12fc 0 (1 - x)) ^ h + h J 12+ko ~ 127 (1 x) ^ 
12 hm 2 ^ + 6 (1 — ug) T b 


(Cla) 

(Clb) 

(Clc) 

(Cld) 

(Cle) 

(Clf) 

(Clg) 

(Clh) 


where 7 , k 0 , S A , S B and T B were given in |B| 


Appendix D: Influence of the cavity size 


In Figs. [ 8 ] and |9j we plot the h dependence of cr c /api ane and A c /Ap; ane , 
f3 = 0.2 and f3 = 0.8, respectively. These plots are qualitatively similar to 
for (3 = 0.5 in Figs. [ 6 ] and [7], and discussed in §VD of the main text. 


for cavity sizes 
those obtained 
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FIG. 8. (a) Critical stress a c normalized by its planar substrate counterpart opi ane , as a function 


of the dimensionless ring thickness h = H/R. Analytical prediction is given by Eq. (21a). (b) 
Critical wavelength A c normalized by its planar substrate counterpart A pi an ei as a function of h. 


Analytical prediction is given by Eq. (21b). Insets are sketches of the wrinkling and global buckling 
modes. The cavity size is /3 = 0.2 and the Poisson’s ratio of the substrate is = 0.5. The legend 
is common to both plots. 
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FIG. 9. (a) Critical stress a c normalized by its plane substrate counterpart api arLe , as a function of 


the dimensionless ring thickness h = H/R. Analytical prediction is given by Eq. (21a). (b) Critical 
wavelength A c normalized by its plane substrate counterpart A pi ane , as a function of h. Analytical 


prediction is given by Eq. (21b). Insets are sketches of the wrinkling and global buckling modes. 
The cavity size is /3 = 0.8 and the Poisson’s ratio of the substrate is us = 0.5. The legend is 
common to both plots. 
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